Schottky-barrier induced spin relaxation in spin injection 
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An ensemble Monte Carlo method is used to study the spin injection through a ferromagnet- 
semiconductor junction where a Schottky barrier is formed. It is shown that the Schottky-barrier- 
induced electric field which is confined in the depletion region and is parallel to the injection direc- 
tion, is very large. This electric field can induce an effective magnetic field due to the Rashba effect 
and cause strong spin relaxation. 

PACS numbers: 72.25Dc, 72.25.Hg, 72.25Rb, 85.75-d 
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Spin injection from a ferromagnetic metal contact into 
a non-magnetic low-dimensional semiconductor struc- 
ture is one of the prerequisites for the realization of 
the next generation high-speed low-power devices based 
on spin degree of freedom. 1 ' 2,3 Notwithstanding the fact 
that many efforts have been devoted to this problem 
experimentally, 4,5 an efficient room-temperature spin in- 
jection is still far away from the horizon. In the mean- 
time, there are many theoretical investigations 7,8,9 on the 
spin injection through the ferromagnet-semiconductor 
junction where a high potential Schottky barrier is 
formed. 6 In these studies the interface is treated through 
various boundary conditions. 7,8 Large (up to 100 %) spin 
injections are reported in these theories. Very recently 
Shen et al. reported a first ensemble Monte Carlo (MC) 
simulation of the spin injection through a Schottky bar- 
rier into a semiconductor quantum well (QW). 9 In this 
study, the Schottky barrier is treated carefully through 
the simulation. Still they reported a substantial spin po- 
larization after the injection to a length scale in the order 
of 1 /im at room temperature without external magnetic 
field. Therefore there must be something missing in the 
theories in dealing with the ferromagnet-semiconductor 
junctions. 

It is noted that the Schottky barrier induces a very 
large electric field parallel to the QW. Such an electric 
field can induce an effective magnetic field due to the 
Rashba effect 10 and can therefore cause a strong reduc- 
tion of the spin polarization after the injection. This ef- 
fect has long been neglected in the literature. A quantita- 
tive estimation of this new relaxation mechanism requires 
an accurate computation of the electric field induced by 
the Schottky barrier which varies strongly with the po- 
sition and is confined in the depletion region. For this 
purpose we apply an ensemble MC simulation to simu- 
late the Schottky barrier and examine the spin relaxation 
induced by this additional relaxation mechanism under 
various conditions. 

We study a ferromagnet-semiconductor diode which is 
one of the elements for many spintronic devices. 11 The 
spin-polarized particles are injected from a bulk ferro- 



magnetic metal into a GaAs QW through a Schottky 
barrier by both thermionic emission and tunnelling in- 
jection, excluding the recombination in the space-charge 
region and the hole injection from the metal to the 
semiconductor. 6 The direction of injection is parallel to 
the QW plane. The electron transport in the QW is 
based on the semiclassical approximation, simply in- 
cluding a "drift" and a "scattering" process: During 
the "drift" process, the spin is influenced by both the 
Rashba 10 and the Dresselhaus 12 spin-orbit interactions. 
The method of the MC simulation has been laid out in 
detail in Ref. 13 for the Schottky barrier simulation, 
in Refs. 14,15 for the spin transport simulation and in 
Ref. 9 for the spin injection simulation. For the inhomo- 
geneous electron distribution in the depletion region, the 
compression/expansion technique presented by Martin et 
al. 16 has to be applied. In this report we do not repeat 
these details except the differences which are addressed 
in the following. 

At finite temperature T, the total current injected from 
a ferromagnetic metal to a semiconductor through the 
Schottky barrier is written as 6 



A*T 



T ms (E x )f m (E)[l-f sc (E)]dE , (1) 



where fee denotes the Boltzmann constant and A* stands 
for the Richardson constant. T ms (E x ) is the tunnelling 
probability through the barrier at the energy E x which 
represents the kinetic energy along the x-direction (the 
injection direction). It is 1 for E x > $b and T ms (E x ) = 
exp{-f J* tp y/2m*[E c (x) - E x ]} for < E x < $ B fol- 
lowing the WKB approximation. <&b represents the 
Schottky barrier height at the metal-semiconductor in- 
terface, xt p is the electron position after tunnelling and 
E c (x) stands for the bottom of the conduction band in 
the semiconductor. f m (E) and f sc (E) are the electron 
distribution functions in the ferromagnetic metal and 
semiconductor separately with E standing for the total 
energy. It is emphasized here that unlike the previous 
works, 9,13 the current and the tunnelling probability are 
only functions of E x , instead of E. After injection, elec- 
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trons start travelling in the QW subject to the spin-orbit 
interactions, the electric field and the electron-phonon 
and possible electron-impurity scattering. The spin-orbit 
interaction is described by H ao (x) — Hr{x) + Hp, with 
the spacial variable-dependent Rashba term 

Hr(x) = j[(<r x ky - <jyk x )£ z - k y (T z £ x (x)] , (2) 



the linear Dresselhaus term = P(k z }(a y k y — a x k x ) 

P{(J x k x k 2 y - 



and the cubic Dresselhaus term H 
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x ). £ x {x) in Eq. (2) is the Schottky-barrier-induced 
electric field (SBIEF). The effective magnetic field in- 
duced by it in Eq. (2) has a significant influence to the 
spin relaxation and has long been overlooked in the lit- 
erature. 

We apply the MC method to study the spin injec- 
tion from the magnetic Fe to (001) GaAs QW through a 
Schottky contact. The well width is 8 nm. The Schottky 
barrier height in the simulation is fixed to be <E>b = 0.72 
eV. 17 We use the following spin-orbit coupling constants: 
(3 = 28 eV-A 3 for the Dresselhaus effect 18 and 7 = 740 
cV-A 2 for the Rashba one. 19 The channel length along the 
spin transport is L x = 2.5 fim. The injection takes place 
at the Fe/GaAs interface at x = 0. As we investigate 
the spin injection from the source, the drain is assumed 
to be in Ohmic contact with the QW. In the figures of 
this paper, we only show the results for the initial 1 fim. 
In the MC simulation, we divide L x into 500 cells and 
choose the time step to be At = 1 fs. To achieve the 
steady transport region, we run the simulation program 
for 10000 time steps and get the results by averaging over 
the last 3000 steps. The initial spin polarization is always 
assumed to be along the x-axis throughout this paper. 
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FIG. 1: Energy of the simulated Schottky barriers for two 
electron densities n at bias V = 0.2 V and T = 300 K. Solid 
curve: n = 10 11 cm~ 2 ; Dotted curve: n — 10 10 cm~ 2 . The 
corresponding electric fields are shown in the inset. 

The simulated Schottky barrier shape, which is de- 
termined by the solution of the Poisson equation and 



the MC simulation of the electron distributions self- 
consistcntly, is shown in Fig. 1 for different electron den- 
sities n in semiconductor QW. An inverse bias voltage 
V = 0.2 V is applied in the simulation, which is in favor of 
the electron injection from the ferromagnetic metal into 
the semiconductors. 6 The large bending near the contact 
indicates the existence of a depletion layer where the elec- 
tron concentration is negligible. It also gives an electric 
field £ x (x) which is shown in the inset of the figure. It 
can be seen that the Schottky barrier becomes thicker 
when the electron density in the QW is decreased. This 
is in consistent with the approximation relation that the 
Schottky barrier width is proportional to n~ 1//2 . 6 Due 
to the change of the shape of the barriers, the SBIEF 
changes also at different electron densities as shown in 
the insect of Fig. 1. 

Because of the large population of the spin-unpolarized 
electrons in the device, especially beyond the depletion 
region, the total spin polarization averaged over all the 
particles at a given position reduces to nearly zero at 
about x = 20 nm. 9,20 We want to get the spin evolution 
of the injected electrons, so our simulation only get the 
spin polarization at each grid averaged over the injected 
spin-polarized electrons. In fact, the spin polarization 
of electrons in the interface of the ferromagnetic metal 
is determined by the spin-dependent density of states of 
electrons in the ferromagnetic contact. Nevertheless, in 
order to investigate spin polarization clearly, we assume 
the injected carrier is S x = 100 % spin polarized first. 

We use |5| = ^JS X + S y + S% to denote the spin polar- 
ization of the injected electrons. Moreover, differing from 
the previous works 9,15 where the electron density is as 
high as 10 12 cm~ 2 , in the present work we only concen- 
trate on the case with density being smaller than 10 11 
cm~ 2 . This is because that when the electron density is 
high, the chemical potential is large compared to fc^T. 
Therefore one should not use the Boltzmann distribution. 
Nevertheless, the MC method treats the scattering semi- 
classically and does not contain any distribution function. 
Consequently, it can only be applied to the problems with 
low electron density. 

In Fig. 2 the spin polarization |5| is plotted as func- 
tion of the position x at temperatures T = 300 and 70 
K without (the red and blue or the dashed and dotted 
curves) and with (the green and pink or the solid and 
chain curves) the SBIEF £ x {x) in the Rashba term [Eq. 
(2)]. The electron density is 10 10 cm~ 2 in the simula- 
tion. We also show the effects of the cubic Dresselhaus 
term to the spin relaxation by performing the simulation 
with (the blue and pink or the dotted and chain curves) 
and without (the red and green or the dashed and solid 
curves) Hp\ It is seen from the figure that the SBIEF 
£ x {x) in Eq. (2) leads to a pronounced spin relaxation 
in the depletion region. The spin polarization is almost 
zero after the depletion layer, in contracts to the previ- 
ous report of a substantial amount at the length scale of 
1 /im. 9 The spin relaxation in all the cases is due to the 
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FIG. 2: (Color online) Spin polarization evolution at different 
temperatures without (the red and blue or the dashed and 
dotted curves) and with (the green and pink or the solid and 
chain curves) SBIEF S x (x) in the Rashba term. The effect of 
the cubic Dresselhaus term is shown by including (the blue 
and pink or the dotted and chain curves) and excluding (the 
red and green or the dashed and solid curves) this term in the 
simulation, n = 10 10 cm -2 . 



D'yakonov and Perel' mechanism. 21 It is further noted 
from the figure that after the fast initial drop of the spin 
polarization in the depletion region, the spin polarization 
also slowly decreases with the position. This is because of 
the spin relaxation induced by the Rasbhba terms from 
the electric field perpendicular to the QW, i.e., £ z , and 
the Dresselhaus terms. This is in contrary to the results 
reported by Shen et al., 9 where they show that the spin 
polarization keeps almost constant beyond the depletion 
region. It is further seen from the figure that the third- 
order Dresselhaus term has marginal effect on the spin 
relaxation, especially at the depletion region where the 
Rashba term is dominant. This is because that the en- 
ergy along the y direction (fcgT) is small, so that fc 2 is 
small compared to fc 2 in the Dresselhaus term. 

We further investigate the effect of the SBIEF to the 
spin injection at higher electron density (but still barely 
in the non-degenerate regime). Curves in Fig. 3 are ex- 
actly corresponding to the cases in Fig. 2 except the elec- 
tron density being 10 11 cm~ 2 , an order of magnitude 
larger. It is seen from the figure that after the deple- 
tion region of Fig. 3, the injected spin polarizations all 
become smaller compared to the corresponding cases in 
Fig. 2. This is due to the enhanced Rashba and Dres- 
selhaus terms at high electron densities. It is seen from 
the figure that the effective magnetic field induced by the 
SBIEF £ x (x) in Eq. (2) again markedly reduces the spin 
injections. It is further noted from the figure that the 

(3) 

cubic Dresselhaus term H y D ' shows larger influence than 
the low electron density case. This is because that the 
SBIEF is much lager at the high electron density case 
(see Fig. 1). This field drives electrons to a much larger 
\k x \ and gives a larger cubic Dresselhaus term. It is also 
seen from Fig. 3 that the effect of cubic Dresselhaus term 
gradually reduces with the decrease of temperature. This 
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FIG. 3: (Color online) Same as Fig. 2 but with the electron 
density n = 10 11 cm -2 . 



is because A: 2 becomes smaller for lower temperature. 
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FIG. 4: (Color online) Comparison of the MC simulation with 
the Blotzmann sampling (a) and the Fermi sampling (b) for 
electrons at relatively high density n = 10 11 cm -2 . T = 300 
K. The meanings of the curves are all the same as those in 
Figs. 2 and 3. 

It is noted that the MC method cannot be applied 
to the strong degenerate (high density) case as reported 
by Shen et al. where the electron density is taken as 
high as 10 12 cm" 2 . 9,15 In the strong degenerate case, the 
electron distribution in the scattering cannot be over- 
looked any more and the MC method fails. Moreover, the 
Boltzmann sampling which is independent of the density, 
should be changed into Fermi sampling. In fact, even for 
the density at 10 11 cm -2 , the non-degenerate approxi- 
mation is already barely valid. In Fig. 4 we show the 
spin injections by using different samplings (Fig. 4(a) for 
Boltzmann sampling and Fig. 4(b) for Fermi sampling) 
at n = 10 11 cm~ 2 , with the scattering still kept to be 
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semiclassical. One can clearly find the marked differ- 
ence. This is because when the Boltzmann sampling is 
used, the energy along the y-axis is fixed in the range of 
ksT, regardless of the density. However, it is much larger 
by the Fermi sampling at high density case. Therefore, 
at high density case, the Rashba and the Dresselhaus 
terms are both weaker from the Boltzmann sampling. 
This leads to a larger spin polarizations. For density at 
10 10 cm -2 reported in Fig. 2, both samplings give the 
same results. 
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FIG. 5: (Color online) Same as Fig. 2 but with the initial spin 
polarization being 40 %. 

In Figs. 2-4, we assume the injected electron spin po- 
larization is 100 %. In fact, the initial spin injection, 
which is determined by the spin-state probability of elec- 
trons in the ferromagnetic contact and can be obtained 
from the microscopic models of ferromagnetic metal, 22 
is about 40 %. Figure 5 shows the results of the spin 
polarization inside the semiconductor QW with electron 



density n — 10 10 cm -2 , when the initial injected electron 
spin polarization is 40 %. The main results are all the 
same as those shown in Fig. 2. As the initial spin polar- 
ization is low, the spin polarization beyond the depletion 
layer is lower than the corresponding case in Fig. 2 and 
reduces to zero more quickly. 

In conclusion, an ensemble MC method is used to simu- 
late the spin-polarized electron injection through a Schot- 
tky barrier and transport in 2D semiconductor QW with 
both Rashba and Dresselhaus spin-orbit coupling. We 
show that the SBIEF not only drives electron to a higher 
momentum states during the injection, which influences 
the spin relaxation via the Dresselhaus and the Rashba 
term, but also provides an effective magnetic field due 
to the Rashba effect. We show that this SBIEF-induced 
effective magnetic field is very strong and gives a pro- 
nounced effect to the spin dephasing at the Schottky 
barrier area. Consequently the spin injection becomes 
almost negligible after the Schottky barrier region. This 
effect has long been overlooked in the literature. More- 
over, this effective magnetic field also provides additional 
relaxation due to the many-body effect 23 which is beyond 
the scope of the MC simulation and will be reported else- 
where. 
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